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Abstract. We calculate the energy spectra and mass compo- 
sition of cosmic rays in energy region (1 -f- 10 8 )GeV/particle 
under the assumption that cosmic rays propagation in the 
Galaxy is described by anomalous diffusion equation. Our 
results and comparisons with experimental data are presented. 



1 Introduction 

The steepening of the all-particle spectrum around 3 • 10 15 eU 
(the "knee") discovered in 1958 (Kulikov and Khristiansen , 
1958) has been the subject of numerous speculations on the 
propagation and acceleration mechanisms of galactic cos- 
mic rays (see, for example, the reviews by Erlykin (1995), 
Kalmykov and Khristiansen (1995), Ptuskin (1997)). The 
usually accepted picture of cosmic ray propagation in the 
interstellar medium is a normal diffusion, which can be de- 
scribed by equation for concentration (Ginzburg and Syrovatskii, 
1964; Berezinsky et al., 1990) 



dN 
~dt 



D(E)AN(r,t,E)+S(r,t,E), 



(1) 



The "knee" is not an intrinsic property of this model under 
the natural physical conditions D(E) ~ E 5 , S(E) ~ E~ p . 
Because of these reasons discussions have mainly been di- 
rected to a search and a justification of breakdown mecha- 
nisms either of the diffusivity D(E) or of the source density 
S(E), However, in spite of considerable theoretical and ex- 
perimental efforts, a model, which can explain 

1 . different values of spectral exponent of protons and other 
nuclei, 



2. mass composition variations at E 
GeV/nucleon, 



10 2 H- 10 5 



was not developed. 

Recently, in our papers (Lagutin et al., 2000, 2001b,c) new 
view of the "knee" problem was presented. It has been shown 
that the "knee" in the primary cosmic rays spectrum is due 
to fractal structure of the interstellar medium, that is another 
regime of particles diffusion in the Galaxy. 

In this paper we consider the propagation of galactic cos- 
mic rays in the fractal interstellar medium taking into account 
that a particle can spend long time in a trap. We demonstrate 
the main results of this new model in a wide energy region 
E ~ (1 -=- 10 8 ) Gev/nucleon. 



2 Model 

Based on the results (Lagutin and Uchaikin, 2001) the cosmic 
ray propagation in the fractal interstellar medium is described 
by fractional diffusion equation. Without energy losses and 
nuclear interactions, the equation for concentration of the 
cosmic rays with energy E generated by sources S(r, t, E) 
has the form 

^ = -D(E, a, 0)D 1 o - (-Ar/ 2 N(r, t, E) + S(r, t, E) 



(2) 



where D(E, a, (3) is the anomalous diffusivity, a and are 
determined by the fractal structure of the medium and by 
the trapping mechanism, correspondingly (see Lagutin and 
Uchaikin (2001)). Dq, denotes the Riemann-Liouville frac- 
tional derivative, (— A) Q / 2 — the fractional Laplacian (Samko 
et al., 1987). 

In the case of point impulse source with inverse power 
spectrum, relating to supernova bursts 



S(r, t, E) = S Q E- p 6(r)0(T - t)Q(t), 0(r) = 



3. the steepening of the all-particle spectrum 



1, T>0, 

0, r<0, 
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the solution of equation (2), found in (Lagutin et al., 2001b) 



is of the form 



N(r,t,E) 



S' E-p 



T -3P/a 



D(E,a,f3) 3 / a 

max[0, t-T] 

x y { t f3) (\r\(D(E,a, ^r^-V^dr, (3) 



where the scaling function ^3 (r), 



* 



(a./?)/ 



W= /4" ) (rT' 3 )gf' 1) (T)T 3 ^ tt rfT, 



(4) 



is determined by three-dimensional spherically-symmetrical 
stable distribution q£*' (r) (a < 2) and one-sided stable dis- 
tribution q[ (t) with characteristic exponent (3 (Zolotarev, 
1983; Uchaikin and Zolotarev, 1999). 

The anomalous diffusivity D(E,a,(3) is determined by 
the constants A and B in the asymptotic behaviour for "Levy 
flights" (A) and "Levy waiting time" (B) distributions: 

D(E, a, j3) oc A(E, a)/B(E, ft). 

Taking into account that both the free path and the probabil- 
ity to stay in trap during the time interval t for particle with 
charge Z and mass number A depend on particle magnetic 
rigidity R, we accept D = (v/c)Do(a, ft)R s . 

Using the representation N — NqE -11 and the property 



To evaluate the parameter ft we have used the results pre- 
sented in the paper (Cadavid et al., 1999), where an anoma- 
lous diffusion of solar magnetic elements have been inves- 
tigated. The authors have shown that the trapping time dis- 
tribution asymtotically takes the form of a Levy distribution 
with spectral exponent ft s=a 0.8. 

Assuming that a trapping mechanism is characterized by 
a kind of self-similarity, one can expect the some value for 
ft in the scales under the consideration . By this reason the 
value (3 = 0.8 is used in our calculations. 

Thus, taking t]e<^e q ~ 2.63 and ??b>b ~ 3.24 we finally 
obtainp«2.9, 5 « 0.27. 

To evaluate the next important parameter — anomalous 
diffusivity Dq(cc, j3), we have used the experimental data on 
the particle anisotropy in the the energy region 10 3 -j- 10 4 
GeV/particle in the framework of the scheme proposed by 
Osborn et al. (1976) and Dorman et al. (1985). For example, 
we find Do « (1 -j- 4) ■ 10~ 3 pc 1 - 7 y^°- s in the case a — 
1.7, (3 = 0.8 for near sources Monogem (r ~ 300 pc, t ~ 
10 5 y), Geminga (r - 300 pc, t ~ 3 ■ 10 5 y), Loop I-IV 
(r ~ (100 -H 200) pc, t ~ (2 -=- 4) • 10 5 y). 

In the model under consideration only one parameter a (1 < 
a < 2) connected with the fractal structure of the interstellar 
medium is found by fit. Extensive calculations of cosmic -ray 
spectra show that the best fit of experimental data may be get 
at a ss 1.7. 



(«,/3) 



jdr = -2 7 rr*^i / 2 ) of the scaling function (Uchaikin 3 Spectra and mass composition 



cm 

and Zolotarev, 1999), one can easy find the spectral exponent 
77 for observed particles: 



The differential flux Ji of the particles of type i due to all 
sources of Galaxy may be separated into two components 



V =P+ ""I 

a 



(5) 



J = J L {r < Ikpc) + J G (r > lkpc). 



(7) 



where 






2-rrr 2 



D(E,a,f3) 2 / a 



.0\-l/a 



j T -^/^^^(\r\(D(E,a,(3y, 

max[0, t-T] ^ 

j ,-3/9/atf («■/») f\r\(D(E, a, (3)TP)- l / a 

max[0, t-T] ^ 



(6) 



Let Eg be a solution of the equation 'E.(E) = 0. One can 
see from (5)-(6) that at E — Eq the spectral exponent for 
observed particles r\ is equal to spectral exponent for particles 
generated by the source: i](Eq) = p. Since the exponent 
Ve<ce = P — S is less than p at E ^C Eo, but the exponent 
Ve»e = P + <V/3 > P at _E > E , E a can be called the 
"knee" energy. 

From experimental values of ?7_e<_Eo an d 7 7e»-E one can 
derive the main parameters of the model (p, S) versus the 
spectral exponent (J3) of "the Levy waiting time": 



ft 



s = (??s»e - '?B«i3o)Y _ ra' p 



VE4CE + S. 



The first component (L) in (7) describes the contribution of 
the nearby sources (at distance r < lkpc) to observed flux 
Ji. The second component (G) is the contribution of the dis- 
tant sources (r > lkpc) to Ji. The similar separation is fre- 
quently used in the studies of cosmic rays (see, for instance, 
Atoyanetal. (1995)). 

The list of nearby sources including 16 supernova rem- 
nants (Nichimura et al., 1979, 1995; Lozinskaya, 1986) is 
used to calculate the L-component. The distant sources are 
supposed to be distributed uniformly both in space and time. 
In this case Jcir > lkpc) ~ E~ p ~ s /P (see Lagutin et al. 
(2001a)). 

Based on this result and (3), we present the differential flux 
in the form: 

J ^) = iJ S N(r i ,t i ,E) + v i C 0i E-P- 6 /P. (8) 

i 
(r<lkpc) 

It is clear from physical point of view that the bulk of 
observed cosmic rays with energy 10 8 -j- 10 10 eU forms by 
numerous distant sources. It means that the observed flux 
at least of protons and He in this energy region must be de- 
scribed by second term in (8). The first term in (8) defines the 
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Fig. 1. Comparison of our calculations of spectra with experimental data: 

Grigirov — Grigorov et al. (1970), 

Ryan — Ryan et al. (1972), 

Minagawa — Minagawa (1981), 

CRN — Grunsfeld et al. (1988), 

HEAO-3 — Engelmann et al. (1990), 

SOKOL — Ivanenko et al. (1990, 1993), 

MSU — Fomin et al. (1991), 

Swordy — Muller et al. (1991) 

Ichimura — Ichimura et al. (1993), 

Zatsepin — Zatsepin et al. (1993), 

Buckley — Buckley et al. (1994), 

Papiny, Lezniak — Wiebel et al. (1994), 

JACEE — Asakimory et al. (1998), 

CAPRICE — Boezio et al. (1999), 

CASA-MIA — Glasmacher et al. (1999), 

AKENO — Youshida et al. (2001), 

RUNJOB — Apanasenko et al. (2001). 



E,Gev 


H 


He 


CNO 


Ne-Si 


Fe 


(In A) 


(A) 


1E+2 


0.53 


0.32 


0.08 


0.05 


0.02 


0.90 


5.36 


3E+2 


0.45 


0.30 


0.12 


0.08 


0.06 


1.20 


8.41 


1E+3 


0.40 


0.29 


0.14 


0.10 


0.08 


1.40 


10.38 


3E+3 


0.36 


0.28 


0.16 


0.11 


0.10 


1.53 


11.65 


1E+4 


0.33 


0.27 


0.17 


0.12 


0.11 


1.65 


12.91 


3E+4 


0.30 


0.26 


0.18 


0.14 


0.13 


1.79 


14.37 


1E+5 


0.26 


0.24 


0.19 


0.15 


0.15 


1.93 


15.95 


3E+5 


0.24 


0.23 


0.20 


0.16 


0.17 


2.06 


17.54 


1E+6 


0.21 


0.21 


0.20 


0.17 


0.20 


2.17 


19.06 


3E+6 


0.19 


0.20 


0.21 


0.18 


0.22 


2.28 


20.44 


1E+7 


0.18 


0.19 


0.21 


0.19 


0.24 


2.37 


21.68 


3E+7 


0.16 


0.18 


0.21 


0.19 


0.26 


2.45 


22.76 


1E+8 


0.15 


0.17 


0.21 


0.20 


0.27 


2.51 


23.70 


1E+8 


0.15 


0.17 


0.21 


0.20 


0.27 


2.52 


23.87 



Table 1. Mass composition of cosmic rays in anomalous diffusion 
model 



spectrum in the high energy region and, as has been shown 
above, provides the "knee". 

We use the spherically symmetric force model of Axford 
and Gleeson (1968) to describe the solar modulation. The 
influence of solar modulation on the particle flux is 



Jmod\-L ) 



T 2 + 2m p c 2 T 



(T + $) 2 + 2to p c 2 (T + $ 



■Jism{T + $), 



where T is the kinetic energy per nucleon, m p is the mass of 
a proton and Jjsm is the interstellar flux (8). The potential 
energy <J>, describing the average energy loss of particle from 
interstellar space to 1 AU, is determined by solar modulation 
parameter <\> : <J> = <j>Z/A. <f> — 750MV is accepted in this 
paper (see Boezio et al. (1999); Menn et al. (2000)). 

The results of our calculation are presented in Fig. 1, and 
Table 1. 



4 Conclusion 

We have considered the propagation of galactic cosmic rays 
in the fractal interstellar medium. The energy spectra of nu- 
clei (H, He, CNO, Ne-Si, Fe) and mass composition have 
been calculated in energy region (1-^10 8 ) GeV/particle under 
the natural physical conditions D(E) - E s , S(E) - E~ p . 
We have shown that the model can explain the different val- 
ues of spectral exponent of protons and other nuclei, mass 
composition variations at E ~ 10 2 -j- 10 5 GeV/nucleon, the 
steepening of the all-particle spectrum. 
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